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Abstract 

We study the quadratic invariants of the elasticity tensor in the framework of its unique 
irreducible decomposition. The key point is that this decomposition generates the direct sum 
reduction of the elasticity tensor space. The corresponding subspaces are completely independent 
and even orthogonal relative to the Euclidean (Frobenius) scalar product. We construct a basis 
set of seven quadratic invariants that emerge in a natural and systematic way. Moreover, the 
completeness of this basis and the independence of the basis tensors follow immediately from 
the direct sum representation of the elasticity tensor space. We define the Cauchy factor of an 
anisotropic material as a dimensionless measure of a closeness to a pure Cauchy material and a 
similar isotropic factor is as a measure for a closeness of an anisotropic material to its isotropic 
prototype. For cubic crystals, these factors are explicitly displayed and cubic crystal average of 
an arbitrary elastic material is derived. 

Key index words: anisotropie elasticity tensor, irreducible decomposition, quadratic invariants 


1 Introduction 

In the linear elasticity theory of anisotropic materials, the relation between the strain tensor Uij and 
the stress tensor the generalized Hooke’s law, is expressed by the use of the elasticity (stiffness) 
tensor 

. ( 1 ) 

In 3-dimensional space, a generic 4th order tensor has 81 independent components. However, due 
to the standard symmetry assumptions for the stress and strain tensors, 

( 2 ) 

the elasticity tensor is left with 21 independent components only. These components are not really 
the intrinsic characteristics of the material because they depend on the choice of the coordinate 
system. Thus, in order to deal with the proper material parameters, one must look for the invariants 
of the elasticity tensor. 


^ijkl _ ^jijkl _ ^klij 


1 
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There are only two linearly independent invariants of the first order of They are usually 

taken as follows, see ISQl, 


Ai = = Cijij and A 2 = j = Cajj . (3) 

Here and subsequently, we use the standard tensor conventions that strictly distinguish between 
covariant and contravariant indices. In these notations, two repeated indices can appear only in up- 
down positions and summation for only two such repeated indices is assumed. The indices of a tensor 
can be raised/lowered by the use of the metric tensor /gki- For instance, the lower components 
of the elasticity tensor are defined as Ciju = gwgjj'gkk'gii'C^ ^ ^ ^ ■ Since in the elasticity literature 
a simplified notation is frequently used, we provide both notations in most cases and relate the 
corresponding quantities by the sign = , as in §• Notice that this shorthand notation is applicable 
only in Euclidean space endowed with rectangular coordinates. 

Quadratic invariants of the elasticity tensor were studied by Ting [25]. He presented two such 
invariants, 

Bi = = CijkiQjki and = C^’^^CBki = CukiCjjki ■ (4) 

Ahmad [T] has contributed the two additional quadratic invariants, 

H 3 = = CukiCjkji and H 4 = = CMuCkjji . (5) 

He also proved that the set of seven quadratic invariants 

1 ^ 1 5 ^ 2 ) ^ 1 ^ 2 , Bi, B 2 , Bs, ( 6 ) 

is linearly independent. Norris |2n) studied the problem of the quadratic invariants and proved that 
the set is complete. It means that every quadratic invariant of the elasticity tensor is a linear 
combination of the seven invariants listed in Q. In order to prove this fact, Norris presented a 
generic quadratic invariant in the form 

T— f, (^ijklf^pqrs /j\ 

^ — J ijklpqrs'^ ^ 1 \> J 

where fijkipqrs is a numerical tensor and provided a detailed analysis of this tensor. It is proven 
that due to the symmetries Q, the most components of fijupqrs vanish and a lot of the remaining 
components are linearly related. This way, exactly seven invariants are left over. In particular, 
Norris demonstrated that a rather natural additional invariant 

Hs = C^^^^Cikji = QjkiCikji (8) 

can be in fact represented as a linear combination of the invariants listed in Q, namely 

B 5 = 2^1 ~ ^ 1^2 + Bi — 2 B 2 + 453 — 2 B 4 . (9) 

In this situation, some principal questions arise: 

• Is there some preferable basis of quadratic invariants? 

• Is there a systematic way to construct such a basis? 
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• Which quadratic invariants can be used as characteristic parameters for different elastic ma¬ 
terials? 

• Which physical interpretation can be given to various quadratic invariants of the elasticity 
tensor? 

In the current paper, we analyze the quadratic invariants problem in the framework of the unique 
irreducible decomposition of the elasticity tensor. In Section 2, we present the principle algebraic 
facts of this decomposition. The key point that such resolution of the elasticity tensor in the 
simple pieces generates the direct sum reduction of the elasticity tensor space. The corresponding 
subspaces are completely independent and even orthogonal relative to the Euclidean (Frobenius) 
scalar product. In this framework, we construct in Section 3 the basis set of seven quadratic 
invariants that emerge in a natural and systematic way. Moreover the completeness of this basis 
and its independence of basis tensors follow immediately from the direct sum representation of the 
elasticity tensor space. We compare this basis to the basis given in (|^. In section 4, we provide 
some applications of quadratic invariants to physics motivated problems. We define Cauchy factor 
for an arbitrary anisotropic material. It can be used as a measure of deviation of a material from 
an analogical Cauchy material. We also prove Fedorov’s relation for an isotropic material closest to 
an anisotropic one. This result follows immediately from the irreducible orthogonal decomposition. 
Correspondingly, we define the isotropic factor that measures the closeness of an anisotropic material 
to its isotropic prototype. In section 5, we study a simplest non-trivial example of cubic crystal. 
It is naturally represented by three independent quadratic invariants. This low-dimensional case 
allows the visualization of the Cauchy and isotropic factors. Corresponding graphs are presented. 
In Conclusion section we present our main results and propose some possible direction for future 
investigations. 

2 Irreducible decomposition 

To describe the irreducible decomposition of the elasticity tensor we first observe two groups acting 
on it simultaneously, see m- For an arbitrary tensor of the range p defined on M”, these are the 
permutation (symmetry) group Sp and the group of rotations 50(n, M). 

2.1 Irreducible decomposition under the permutation group 

The symmetry group S '4 provides permutations of the elasticity tensor indices. The decomposition 
of under this group is described by two Young diagrams: 



( 10 ) 


All other 4-the order Young’s diagrams are non-relevant in our case due to the original symmetries 
([^ of , The left-hand side of (10) represents a generic 4th rank tensor. On the right-hand 
side, two diagrams describe two tensors of different symmetries. The explicit expression of these 
tensors can be computed by the use of the corresponding symmetrization and antisymmetrization 
operations. For the elasticity tensor, the result is obvious. The first (row) diagram represents the 
totally symmetric tensor. The second (square) diagram represents an additional tensor that can be 
considered merely as a remainder. 
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As a result, we arrive at the unique and irreducible decomposition of the elasticity tensor under 
the action of the group ^ 4 : 

(jijkl _ gijkl _|_ ^ 

with 

gijkl fjiijkl) _ 1 ^(jijkl _|_ Qiklj _|_ fjiljk'^ ^^2) 

and 

j^ijkl Qijkl _ (j(ijkl) _ \ ^Qijkl _ ^ilkj _ fjiklj'^ 

We use the standard normalized Bach parentheses for symmetrization and antisymmetrization of 
indices. We observe that every antisymmetrization of gives zero and every symmetrization 
preserves this tensor. As for the second part its total symmetrization vanishes. Con¬ 

sequently, we have a useful symmetry relation 

j^ijkl _|_ j^iklj _|_ jjiljk _ q ^ 

As it is shown in m and m, the equation A^^^^ = 0 describes the well-known Cauchy relation. 
Thus, we call the Cauchy part and A^l^^ the non-Cauchy part of the elasticity tensor . 
Observe the main algebraic properties of this decomposition: 

• The partial tensors and satisfy the minor symmetries, 

S[ij]kl ^ gijlkl] ^ Q J^[ij]kl ^ J^ij[kl] ^ Q ^ 

and the major symmetry, 

gijkl^gklij j^ijkl ^ ^klij _ 

Thus, these partial tensors can themselves serve as elasticities of some hypothetic material. 

• Moreover, any additional symmetrization or antisymmetrization preserves the tensors 
and or nullifies them. 

• The decomposition © is preserved under arbitrary linear transformations. Thus, it can be 
referred to as irreducible GL(3, M)-decomposition. 

• The irreducible decomposition of the tensor C provides the decomposition of the corresponding 
tensor space C into a direct sum of two subspaces § C C (for the tensor S) and A C C (for the 
tensor A), 

e = S©A. (17) 

In particular, we have 

dimC = 21, dimS = 15, dimA = 6. (18) 

• The irreducible pieces 5®-^^^ and A®-^^^ are orthogonal to one another in the following sense: 

= ( 19 ) 

Indeed, 


( 20 ) 
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• “Pythagorean theorem:” The Euclidean (Frobenius) squares of the tensors 


/o 2 /^ijkl 

^ = 0 Cijkl , 

02 oiikl 0 

.b = bijki , 

^2 ^ A^jki^,,^^ 

(21) 

satisfy the relation 

= S^ + A^. 


(22) 


2.2 Irreducible decomposition under the rotation group 

We are looking now for the decomposition of the elasticity tensor under the action of the group of 
rotations. The following fact is due to the classical theory of invariants: Relative to the subgroup 
SO{3) of G'L(3), a basis of an arbitrary system of tensors coincides with a basis of the same 
system with the metric tensor added m- We use the Euclidean metric tensor gij. In rectangular 
coordinates, gij = diag(l, 1,1). 

We start with the totally symmetric Cauchy part From the contraction of with the 

metric tensor, we construct a unique symmetric second-rank tensor 

5'^' := ^ (Qjkk + 2Cikkj) (23) 

and a unique scalar 

S := gijS^^ = Siikk = - {Ciikk + ‘^Cikki) ■ (24) 

We denote the traceless part of the tensor as 

:=S^^ -^Sg^\ with gijP^^=0. (25) 

Now we turn to the decomposition of the tensor We denote the two subtensors 

{i)gijkl _ ppiijgkl) . ( 26 ) 


It can be checked now by the straightforward calculations that the remainder 

_ i^ijkl _ [l)i^ijkl _ {2)i^ijkl 

is totally traceless if and only if 


1 

Oi —, 

5 




(27) 


(28) 


Hence, we obtain the decomposition of the totally symmetric tensor into the sum of three 
independent pieces: 

gijkl _ {l)gijkl _|_ {2)gijkl _|_ {3)gijkl (2g^ 


where 


{l)gijkl ^ = ^5 (gijgkl ^ gikgjl ^ gilgjk\ _ (3Q^ 

5 15 V / 

{2)gijkl _ ^piijgkl) _ 1 ^pijgkl _|_ pikgjl _|_ pilgjk _|_ pjkgil _|_ pjlgik _|_ pklgij^ _ ( 3 ^^ 


{3)pijkl _ pijkl 


and 


( 32 ) 
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These pieces are unique and invariant under the action of the group SO{3). Moreover, the cor¬ 
responding subspaces and are mutually orthogonal. Indeed, for A,B = 1,2,3 with 

A^B, 

= 0 . (33) 

This fact follows immediately from the tracelessness of the tensors and 

Consequently, the vector space § of the totally symmetric tensor is decomposed into the 
direct sum of three subspaces 

§ ^ (I)g 0 (2)g 0 (3)g (34^ 

with the corresponding dimensions 

15 = 1©5©9. (35) 

We turn now to the second part of the elasticity tensor. The irreducible piece is a fourth 
rank tensor with 6 independent components. It is quite naturally that it can be represented as a 
symmetric second-rank tensor, see mm, [T3] . and [H] for detailed discussions. We dehne 

^ran ■— '^^mil^njk-^ ^ , (36) 

where Cijk = 0, ±1 denotes the 3-dimensional Levi-Civita permutation pseudo-tensor. Consequently, 
Amn is a symmetric tensor, Amn = Anm, and we have 

^ijki ^ _ ( 37 ^ 


For a proof of this proposition it is enough to substitute (36) into (37) and to apply the standard 
relations for the Levi-Civita pseudo-tensor. 

In order to decompose the non-Cauchy part A^l^\ it is convenient to use its representation by 
the tensor density Aij. We denote 

A := g^^A^n • (38) 


By using the relation 


9 ^mil^njk 9ij9lk 9ik9jl 


(39) 


and Eq.(37), we derive 


^ — 3 {gijOlk 9ik9jl) ^ — 3 i-^iikk -^ikik) — 3 {Ciikk Cikik) ■ 

The tensor density Aij can be decomposed into the scalar and traceless pieces: 


(40) 


Aij — Qij -\- ^Agij , 


(41) 


where the traceless piece is given by 


Qij :— Aij ^Agij 


(42) 


Substituting (41) into (37), we obtain 


J^ijkl _ (l)jjijkl _|_ (2)j^ijkl 


( 43 ) 
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where the scalar part is given by 

(2g^lg>^^ - gilgjk - g^^gl>^^ (44) 

and the remainder reads 

We recall that the tensor Qmn is symmetric and traceless. Since the product of two Levi-Civita 
pseudo-tensors is represented by the determinant of the metric tensor the latter equation can 
be rewritten as 

{2)j^ijkl gikQjl ^ gjkQil ^ gilQjk ^ gjlQik _ 2gklQij _ 2gilQf^l . (46) 

The decomposition given in Eq.(|43[) is unique, invariant, and irreducible under the action of the 
rotation group 50(3, M) and of the permutation group S 4 . 


Correspondingly, the vector space A of the tensor is irreducibly decomposed into the direct 
sum of two subspaces 



yi = (byi © (2)yi ^ 

(47) 

with the corresponding dimensions 

6 = 1©5. 

(48) 

Since the trace of equals zero, 

these subspace are orthogonal to one another, 




(49) 


Collecting our results, we formulate the following 


Theorem 1. Under the simultaneous action of the groups and 50(3, M), the elasticity tensor is 
uniquely irreducibly decomposed into the sum of five parts 


5 

^ijkl = y. {A)fjijkl 


(50) 


This decomposition corresponds to the direct sum decomposition of the vector space of the elasticity 
tensor into five subspaces 

e = © (5)g^ ^ (54^ 

with the dimensions 

21 = (1 © 5 © 9) © (1 © 5) , (52) 

The irreducible pieces are orthogonal to one another: For Afi^B 

= (53) 

The Euclidean squares, = CijkiC^^^^ and with A = 1, • • • ,5, fulfill the 

“Pythagorean theorem:” 

C2 ^ (2)^2 (3)(^2^ J'(4)©,2 (5)©,2^ _ 


(54) 
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2.3 Irreducible decompositions 


The decomposition (51) involves two scalars S and A, two second order traceless tensors Pij and Qij, 
and a fourth order totally traceless tensor Rijki- Exactly the same types of tensors emerge in the 
harmonic decomposition that is widely used in elasticity theory. Such decomposition is generated 
by expressing the partial tensors in term of the harmonic polynomials, i.e., the polynomial solutions 
of the Laplace equation. The corresponding tensors are required to be completely symmetric and 
totally traceless. As it was demonstrated by Backus [3], such a harmonic decomposition is not 
applicable in general in a space of the dimension greater than three. 

The most compact expression of this type was proposed by Cowin [8j , 

C^ikl ^ ^gijgkl ^ ^ ^gikgjl ^ gilgjk^ ^ ^ ^ ^ ^ ^ ^ ^ ^ i ^ + 

gil^jk gjk^il gik^jl'^ ^ijkl _ 

An alternative expression was proposed by Backus [3]. It reads, see [1], 

(Jijkl ^ fjijkl ^ (^H^jgkl ^ fjikgjl ^ ^ilgjk ^ ^jkgil ^ ^ + 


(55) 




il jk 




+ 




(56) 


Let us compare these two expressions. First we observe that two totally traceless tensors must be 
equal to one another, As for the scalar terms in Eq.(55), they are merely linear 

combinations of the corresponding terms in Eq.(56). Indeed, it is enough to take 


a = H + h and b = H — {l/2)h. 


(57) 


Quite similarly, the traceless second order tensors of Eq. (55) are linear combinations of the corre¬ 
sponding terms of Eq.(56) with the identities 

^ki ^ jjki ^ j^ki ^ji ^ jjji _ ^i/2)hA . (58) 

Both decompositions are irreducible under the action of the rotation group. The key difference 
between the two is that that the decomposition of Backus is also irreducible under the action of 
permutation group. Cowin’s decomposition is reducible in this sense. 

Let us compare now the harmonic decomposition Eq.(56) to our decomposition as it is given in 
Eq.(50). We immediately identify 

{i)Sijki ^ + 

(2)gijkl 
(i)gijkl 

{l)j^ijkl 


[2)j^ijkl _ 


Riig^l + W^gA + gil ^ jjjlgik ^ ^klgij ^ 

jjijkl 

h - \g^’^g^^ 

f^ijgkl ^ f^klgij _ \h3lgik _ Ih^kgjl _ lykgil _ ^hUglk _ 


(59) 

(60) 
(61) 

(62) 

(63) 
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We can straightforwardly derive the relations 


"=15®- 


— _ptj 

7 


j^ijkl _ pijkl 


and 




= -2Q 




(64) 


(65) 


Thus the two decompositions are equivalent. The difference is that in Eq.(50) the partial tensors 


are identified as elasticities themselves. These partial tensors generate the minimal direct sum 
decomposition of the elasticity tensor space. Moreover, the corresponding subspaces are mutually 
orthogonal and the squares of the tensors satisfy the “Pythagorean theorem”. We will see in the 
following section how those properties can be applied to the problem of quadratic invariants. 


3 Linear and quadratic invariants 

3.1 Linear invariants 

Each linear invariant of the elasticity tensor can be represented as 

I = , 


( 66 ) 


where hijki is a tensor. Using the decomposition (51), it can be rewritten as a sum of irreducible 
parts with different leading coefficients 




(67) 


1=1 


We observe that the tensor hijki can be constructed only as a product of two components of the 
metric tensor, hij^i ~ {g ■ g)ijki- Since the tensors Pij, Qij, and Rijki are totally traceless they do 
not contribute to the sum in Eq. (67). Consequently we are left with 


I = aS + l3A. (68) 

Hence, every linear invariant is represented as a linear combination of the two basic linear invariants 

(69) 


Li - S - g {Cukk + 2Cikik) , ^2 - A - ^ {Cukk C,kik) . 


Their independence can be seen from their explicit expressions. It follows also from the fact that S 
and A are related to two different subspaces and We readily obtain the expression of the 
linear invariants given in Eq.([^, 


Cikik — S — A, Ciikk — S + 2A . 


( 70 ) 
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3.2 Quadratic invariants 

Norris |2U] presented a generic quadratic invariant of the elasticity tensor in the form 


f.... (~<ijkl(-,pqr 
J ijklpqrs^ ^ 


(71) 


Here fijkipqrs is a numerical tensor. Using the irreducible decomposition (51) it can be rewritten as 


5 


/ ^ JzjkLpqrs ^ ^ 

1J=1 


(72) 


The tensors with the components fijkipqrs can be constructed only as a product of four compo¬ 
nents of the metric tensor, fijkipqrs ~ {d ' 9'9 ' g)ijkipqrs- Using the traceless property of the 
tensors Pjj, Qij, and Rijki, we can show that the quadratic invariants can be chosen uniquely as 



Zi = S\ 

Z2 = AS, Z3 = A\ 


(73) 

Z 4 = Pijp'^ , 

Z 5 = RjQ '^, 

Zq = QijQ^^ , and 

Z 7 = RijkiR^^'^'. 

(74) 


It is clear that this set of invariants is complete. Indeed, when two tensors in (72) are irreducibly 
decomposed in the form (51) with an arbitrary tensor fijkipqrs constructed from the metric tensor 
and numbers, only the terms (73, 74) can appear. Moreover, these invariants are independent. It is 
due to the fact that they are taken from independent and even orthogonal subspaces of the elasticity 
tensor space. 


3.3 Relations between two sets of qnadratic invariants 

Let us display the quadratic invariants of the set ([^ in terms of Zj. We present the details of 
calculations in the Appendix. For the quadratic invariants constructed from the linear ones, we 
have straightforwardly 


Al = {S-Af = Zi-2Z2 + Z3, (75) 

A 1 A 2 = (S-A)(S + 2A) = Zi + Z2-2Z3, (76) 

Al = (5 + 2A)2 = + 4 Z 2 + 4^3 ■ (77) 

For the first invariant of Ting, we write 

5 5 

Bi = C^^’^^'Cijki = E {I)(-iijkl E ^^^C^jki . (78) 

i=i j=i 

Due to the orthogonality of the set, we are left here with 

P, = (1)^2 ^ (2)^2 ^ (3)^2 ^ (4)^2 ^ (5)^2 (79) 


Let us list the expressions of these invariants 

(1)^2 ^1^2, 

5 


( 3)52 ^ Rijki^ 


( 2)^2 ^ ^pijp^. , 


(80) 
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and 

{l) A'i _ a2 ( 2 ) 42 _ 

Consequently, 

1 6 

Bi = -Z\ + Zs + -^4 + 12^6 + -^7 • 

5 7 

The second invariant of Ting takes the form 

B2 = -Zi + -Z2 + -Z^ + Zi — 4Z5 + dZg ■ 
The first and the second invariants of Ahmad read 

-B3 = + 2-^2 — + -^4 — -^5 — 2^6 


(81) 

(82) 


(83) 


(84) 


and 

-B 4 = -Zi — -^2 + 2^3 + -^4 + 2^5 + Zg ) 


(85) 


respectively. From these expressions we see that the set of invariants is complete and the 
invariants are independent. The same is true for our set Zi. We calculate also the additional 
invariant of Norris: 


55 = C^^^'Cikji = ^Zi - ^Z3 + ^Z4 - 6Ze + ^7 . 


( 86 ) 


Substituting the expressions (73 74) we obtain the formula ([^ of Norris. 


4 Applications: Invariants as characteristics of materials 

Since invariants of the elasticity tensor are independent of the coordinate system used in specific 
measurements, they can be used as intrinsic characteristics of the materials. It is clear that linear 
independence is not enough for this goal. Indeed, although the invariants Zi, Z 2 ) and Z 3 are linear 
independent, they are related by a quadratic relation Z| = Z 1 Z 3 . We will show that an intrinsic 
meaning can be assigned to the five invariants that correspond to different direct subspaces of the 
elasticity tensor space, namely 

{Zi , Z 3 , Z 4 , Zq , Zyj . (87) 

All these invariants are positive. 

4.1 Cauchy relations and Cauchy factor 

In the early days of the elasticity theory, Cauchy formulated a molecular model for elastic bodies, 
based on 15 independent elasticity constants. In this way 6 constraints, called Cauchy relations 
were assumed. A lattice-theoretical analysis shows, see m, m , that the Cauchy relations are 
valid provided the following conditions hold: 

• The interaction forces between the molecules of a crystal are central forces; 

• each molecule is a center of symmetry; 

• the interaction forces between the building blocks of a crystal can be well approximated by a 
harmonic potential. 
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More recent discussions of the Cauchy relations can be found, e.g., in [1], [3], [4], or [7]. Different 
compact expressions of the Cauchy relations can be found in literature. For instance in m. they 
are presented as 

^iijk _ (Jijik ^ Q _ 

An alternative form is widely used, see m. m. uni, 0 , 

^ijkl _ (Jikjl ^ Q _ (gg^ 

The irreducible decomposition technique m yields 

= 0, or Qij = 0. (90) 


As it was demonstrated experimentally already by Voigt, the Cauchy relations do not hold even 
approximately. Thus, elastic properties of the generic anisotropic material is described by the whole 
set of 21 independent component. In fact, the situation with the Cauchy relations is much more 
interesting, see [12] . One can look for the deviation of the elasticity tensor from its Cauchy part. 
As it was pointed out by Haussiihl this deviation, even being a macroscopic characteristic, 
can provide some important information about microscopic structure of the material. To have 
such deviation term we must have a unique proper decomposition of the elasticity tensor into two 
independent parts that can be referred to as Cauchy and non-Cauchy parts. 

In [T2|, the deviation from the Cauchy part was presented by the value of the corresponding 
combination given in the left hand side of Eq.(88). Such way of expression is valid only in the 
case when the non-Cauchy part is presented by only one component. Moreover this expression is 
dimension-full and depends on the choice of the coordinate system. 

With the use of the quadratic invariants we can introduce an invariant characteristic of deviation 
of a material from its Cauchy prototype. Due to the “Pythagorean theorem” (51), we define the 
dimensionless quantity, which we will call the Cauchy factor 


3”Cauchy — 




(91) 


Evidently, 0 < Tcauchy < 1. A pure Cauchy material is determined by Tcauchy = 1- For Tcauchy = 0, 
we have a hypothetic material without Cauchy part at all. Comparing two materials, we must 
conclude that a material with higher Cauchy factor has a microscopic structure closer to spherical 
symmetry. 


4.2 Fedorov’s problem 

In linear elasticity for anisotropic materials one must deal with a big set of elasticity constants. But 
in some problems, the elastic body can only be slightly different from an isotropic one. Fedorov 
in a classical book on the propagation of elastic waves in anisotropic crystals, has demonstrated 
how the anisotropic elastic tensor can be averaged over the 3-dimensional spatial directions in order 
to hnd some kind of isotropic approximation. 

Recently Norris |21) took up this program and defined an Euclidean distance function for solving 
the Eedorov problem in a novel way. He succeeded in doing so and even extended the formalism for 
averaging the given set of elastic parameters relative to less symmetric classes. The corresponding 
procedure can be outlined as follows: 
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• For a given 4th order elasticity tensor, one constructs the corresponding 2nd order Christoffel 
tensor, which is quadratic in the wave vector. 

• One consider the norm of the difference between the given anisotropic Christoffel tensor 
and a generic isotropic one. The isotropic Christoffel tensor is taken with two unknown 
parameters. 

• Since Christoffel tensor is quadratic in the wave vector n, the mentioned norm is quartic n 
likewise. Its average is computed in space directions. 

• The resulting expression is left to be a function of two isotropic parameters. Its minimization 
is applied and the resulting pair of isotropic parameters is derived. 


The result of this consideration is given in |21j as 

K = g Ciijj, ^ ~ IQ ~ ^ ■ (92) 

In 120], [21] and (TH], Norris explained that minimizing a Euclidean distance function is equivalent 
to projecting the tensor of elastic stiffness onto the appropriate symmetry. 

Let us consider an elasticity tensor of 21 independent components. It is irreducibly decomposed 
to the sum of five independent pieces. Two scalar pieces, namely a,nd has a special 

property: They are invariant under arbitrary SO{3) transformation. Using the direct sum of the 
corresponding subspaces, we can construct a subspace 

ISO = (93) 


This 2-dimensional subspace is SO(3) invariant and orthogonal to all other subspaces of an arbi¬ 
trary elasticity tensor. Evidently it must be identified as an isotropic part of an elasticity tensor. 
Consequently we constructed an isotropic part of the generic elasticity tensor in the form 


(iso 


)(jijkl 


(l)gijkl (l)jyijkl 


(94) 


or, explicitly, 

iiso)cijkl ^ (^gijgkl ^ gikgjl ^ ^ ^ ^ (^^g^jgkl _ 

Let us compare this expression to the standard representation of an isotropic material in terms of 
the Lame moduli A and p 

{iso ) ^ijkl ^ ^gijgkl ^ ^ ^ 

We derive the effective Lame moduli for an anisotropic material 

S — “ {Ciikk T , A — — {Ciikk . 


Recall that 


(98) 
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Substituting these expressions into equations (97), we derive 




IJIJ ? 


* 1 1 

~ IQ ^ Ciijj 


(99) 


With the bulk constant k = X + (2/3)/r we recover both expressions given in Eq.(92). 

In order to express the deviation of the given anisotropic material from its effective isotropic 
prototype, one uses the distance between two tensors. This quantity is dimensionful and depends 
on the average magnitude of the elasticity tensor. Instead, we define the isotropy factor of an 
anisotropic material in the form 


Tisn = 


(iso) (-72 




In terms of the constants S and A it reads 




Tiso = 


52 + 5^2 


Z\ + 5 Z 3 


y Zi + 5^3 + (30/7)y4 + 60Zq + 52^7 


( 100 ) 


( 101 ) 


We observe that 0 < Tiso <1- It is equal to one for pure isotropic materials and equal to zero for 
some hypothetic material without isotropic part, i.e., in the case when the effective Lame moduli 
vanish. 


4.3 Irreducibility factors 

As a natural extension of the Cauchy and the isotropy factors described above, we introduce dimen¬ 
sionless numerical factors that describe the contribution of the irreducible pieces to the elasticity 
tensor. For the 5 irreducible parts with / = 1, • • • , 5, we define the irreducibility factors 


kJ irr 




( 102 ) 


W Aijkl^^) Adkl + 


In particular, the Cauchy factor is expressed as 

Cauchy 

and the isotropy factor as 
3^iso — \ 


ijkr 








ijkr 






(103) 


(104) 
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5 Cubic crystals 

5.1 Definition 

Cubic crystals are described by three independent elasticity constants. In a properly chosen coor¬ 
dinate system, they can be put, see Nayfeh m, into the following Voigt matrix: 


(105) 


-(-nil 

^1122 

^1133 

(-1123 

^1131 

(-1112- 



^12 

(-12 

0 

0 

0 ■ 

* 

^2222 

^2233 

(-2223 

(-2231 

(-2212 


* 


(-12 

0 

0 

0 

* 

* 

^3333 

(-3323 

(-3331 

(-3312 


* 

* 


0 

0 

0 

* 

* 

* 

(-2323 

(-2331 

(-2312 

— 

* 

* 


(-®® 

0 

0 

if: 

* 

* 

* 

(-3131 

(-3112 


* 

* 

* 

* 

c®® 

0 

if: 


* 

* 

* 

(-1212 


* 


* 

* 

* 

c®® 


Taking into account the multiplicities of the elements exhibited in ( |105[ ), we calculate 

= 3 + 6 + 12 . 

5.2 ^4 -decomposition 


We decompose (105) irreducibly and find the Cauchy part 


gijkl _ 


-51111 

^1122 

^1133 

51123 

51131 

511121 


a 

/3 

/? 

0 

0 

O' 

* 

^2222 

^2233 

52223 

52231 

52212 



a 

/3 

0 

0 

0 

* 

* 

^3333 

53323 

53331 

53312 



* 

a 

0 

0 

0 

* 

* 

* 

52323 

52331 

52312 


* 

* 

* 

/? 

0 

0 

* 

* 

* 

* 

53131 

53112 


* 

* 

* 

* 

/3 

0 

* 

* 

* 

* 

* 

51212 


* 

* 

if: 

* 

* 

/3. 


where 


a = /3 = J + 2C®®) , 

o 


66^2 


The square of this tensor takes the value 

5^ = + 18/3^ = 3 + 2 + 2C®®) 

The non-Cauchy part is represented by 




-^1111 

^1122 

^1133 

^1123 

^1131 

^11121 


0 

27 

27 

0 

0 

0 ■ 

* 

^2222 

^2233 

^2223 

^2231 

^2212 


* 

0 

27 

0 

0 

0 

* 

* 

^3333 

^3323 

^3331 

^3312 


* 

* 

0 

0 

0 

0 


* 


^2323 

^2331 

^2312 


* 


* 

-7 

0 

0 


* 

* 

=t= 

^3131 

^3112 


* 

* 

* 

* 

-7 

0 

* 

* 

* 

=t= 

=t= 

^1212 


* 

* 

* 

* 

=1 

-7. 


where 


Its square reads 


7 = 3 - C®®) 


= 367^ = 4 - C®®)^ . 


(106) 


(107) 


(108) 

(109) 


( 110 ) 


( 111 ) 

( 112 ) 
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We can straightforwardly check the identities 


^•5 = 0, 


and 

= S'^ + A^ = + 18^2 ^ 35^2 _ 

The Cauchy factor is expressed as 

^ _ I q;2 + 6^02 

TCauchy - y + 6/32 + 1272 ’ 


or 


3”Cauchy 


For the Cauchy relation we have 


I 3(C11)V2(C12 + 2C66)2 
3(C11)V6(C12)2 + 12(C66)2 ■ 


7 = 0 , or . 


(113) 

(114) 

(115) 

(116) 

(117) 


In this case, Eq.(115) yields Tcauchy = 1- 

For cubic crystals, the Cauchy factor in Eqs.( 115|116 ) depends on two independent parameters 
and thus allows 3-dimensional visualization, see Fig. 1 and Fig 2. 



Figure 1: Functional dependence of the Cauchy factor Tcauchy with respect to the parameters (3/a 
and 7 /a. The pure Cauchy materials are depicted by the straight line lying on the axis 7 = 0. 


5.3 S'(9(3)-decomposition 


Let us determine the SO{3) components of the totally symmetric tensor. Using Eq.(23) we derive 
from (107) 


Consequently, 


Sij 


= (a + 2l3)g^^ = 
S = 3(a + 2/3) 


^ 2^12 ^ 4 ^ 66 j gij ^ 

(118) 

3C^^ + 2C^2 ^ 4,766 ^ 

(119) 

















y. Itin 


Quadratic invariants of the elasticity tensor 


17 




Figure 2: Functional dependence of the Cauchy factor 3“cauchy with respect to the parameters 
^ 12/^11 The pure Cauchy materials are depicted by the straight line = C®®. 


and 

= 0. 

The total traceless remainder can be expressed as 


j^ijkl _ 


The two non-zero quadratic invariants read 

= 9{a + 2/3)2 = (3C^^ + 

and 


We check the formula 


i?2 = - 3/3)2 ^ I /^ll _ ^12 _ ^66^2 

5 5 ^ ^ 


1 


52 = -S^ + R^, 
5 


( 120 ) 



^1122 

^1133 

^1123 

^1131 



[2 

-1 

-1 

0 

0 

0 ] 

* 

j^2222 

^2233 

^2223 

i?2231 

^2212 



2 

-1 

0 

0 

0 


* 

^3333 

j^3323 

^3331 

^3312 

a — 3/3 

* 

* 

2 

0 

0 

0 

* 

* 

* 

j^2323 

^2331 

^2312 

5 

* 

* 

* 

-1 

0 

0 


* 

* 

* 

^3131 

^3112 


* 

* 

* 


-1 

0 

* 


* 

* 

* 

^1212 



* 

* 

* 


- 1 _ 


( 121 ) 

( 122 ) 

(123) 

(124) 


that must hold with the correspondence with the expressions ( 80|81 ) of the first Ting’s invariant. 
Let us turn now to the asymmetric part Eq.(llO) yields 


Consequently, 


A*-’' = 275 *^'. 

A = 67 , = 0 . 


(125) 

(126) 


Thus, 


= A ^ = 3672 = 4 (Ci 2 _ (766^2 ^ 


( 127 ) 
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Accordingly, the isotropic factor takes the form 


T- — 

isn — 


(l/5)5'2 + A2 /(3/5)(a +2,5)2+ 1272 




a2 + 6/32 _|_ 1272 


(128) 


or 


Ifian = 


\ 


(3C11 + 2C12 + 4C766)2 + 20 


15 


+ 2(C12)2 +4(C66) = 


(129) 


It is well known that an isotropic medium can be described as a cubic crystal with C®® = (1/2) (C^^ — 
( 7 ^ 2 ), It ig equivalent to a = 3/3. When this relation is substituted into (128), we obtain Ibiso = 1- 


In Fig. 3 and Fig 4. we present the functional dependence of the isotropy factor for a cubic crystal 
with respect to the homogeneous fractions of its parameters. 



Figure 3: Functional dependence of the isotropy factor IFiso is depicted with respect to the variables 
Y/a and y/a. The straight line with fd/a = 1/3 corresponds to the pure isotropic materials. 



0 


| 5 /a 



Figure 4: Functional dependence of the isotropy factor Tiso with respect to the variables 
and jC^^. The straight line C®® = (1/2)(C'^^ — C^^) correspond to the isotropic medium. 







































y. Itin 


Quadratic invariants of the elasticity tensor 


19 


5.4 Cubic averaging 


We consider now a problem of averaging of an arbitrary elastic material by a cubic symmetry pro¬ 
totype. Let in some coordinate system a material be described by a full set of 21 elastic parameters 
(jijki ^ "We are looking for a cubic crystal that is mostly close to our material with respect to the 
Euclidean metric. We assume that in the same coordinate system the cubic crystal is represented 
in the canonical form (105). Since for cubic crystal = 0, its elasticity tensor is left 


where 


II 

0 “ 

(1) qijkl 1 (3) qijkl , (1) Aijkl 
•^cub “^cub Wub’ 

(130) 

(1) qijkl 
‘^cub ~ 


(131) 


( 25 V'- 5V"- <7* V) , 

(132) 

(3) gijkl _ 


(133) 


Here is a set of parameters (not a tensor) that is represented in the 6x6 notations by the 
constant matrix in the right hand side of Eq.( 121| ). Thus m,k,n are our unknown variables. The 
square distance between two elasticity tensors is given by 


-c: 


'cuh) =Ki{S- mf +K 2 {A- nf + + P 3 , 


(134) 


where Ki, K 2 , are positive numerical constants. This expression reaches its minimal value for 


m = S, n = A, 


and 

We have 




k = 




Using the expression of listed in Eq.(121), we derive 

_)_ j :^2222 _|_ ^3333 


k = 


6 


Consequently the parameters of the cubic crystal prototype take the values 

a = -5 + 2fe = - (C11 ^i + C2222 + c3333) _ A5 

5 3 15 

= u, 

6 


(135) 

(136) 

(137) 

(138) 

(139) 

(140) 

(141) 


7 
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6 Conclusion 

In the framework of the irreducible decomposition of elasticity tensor, we studied the problem of 
its quadratic invariants. Since this decomposition is orthogonal, the invariants emerge in a natural 
and systematic way. Their independence and completeness follow straightforwardly from the direct 
sum decomposition of the tensor space. For arbitrary anisotropic materials, we defined the Cauchy 
factor as a dimensionless measure of a closeness to a pure Cauchy material. Quite similarly, we 
defined isotropy factor as a measure for a closeness to an isotropic prototype of a given material. 
The irreducible factors are defined in order to characterize the contributions of different irreducible 
parts of an anisotropic elasticity tensor. This formalism can be useful for various elasticity problems: 

• Elasticity wave propagation [2]; 

• complete set of anisotropy invariants, see [28] : 

• material symmetries and wavefront symmetries 0, 0; 

• averaging of anisotropic material by a higher symmetric prototype [18], [21] . 
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A Calculating the relations between invariants 

A.l The first Ting invariant 

For the first invariant of Ting, we write 

5 5 

E 

I=l J=1 

Due to the orthogonality of the decomposition, we are left with 

+ (2)^2 ^ (3)^2 ^ (1)^2 ^ (2)^2 

We calculate step by step 

(2) ^2 _ _j_ pilgjk _j_ pjkgil _|_ pjlgik _|_ pklgij^ 

(3) q2 jyijkl T} _ <7 

J ^ ^ijkl — ^7 , 


and 


(142) 

(143) 

(144) 

= ®P‘^R., = ®y4,(145) 

(146) 

(147) 

(148) 


(1) yl2 ^ J_^2 / 2gdgkl _ gilgjk _ gikgjl\ ^ = ^2 ^ 

36 V / 

( 2 ) ^2 ^ 1 ^ ^ ^ = 12^6 • ( 149 ) 
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Consequently, the first invariant of Ting reads 



Bl — -Z\ + ^3 + -^4 + 12^6 + . 

5 7 

(150) 

A.2 

The second Ting invariant 


For the second invariant of Ting, B 2 , we first calculate the trace 



C\ki — - 9 ^^ {^Sijki + . 

(151) 

Here, 

s'’Smi = Sa = Pa + \sg9i, 

(152) 


9^^ — 9^ {,9ikQjl H“ 9jkQil “1“ 9ilQjk 9jlQik ‘^9klQij ‘^9ijQkl) — ^Qkl ? 

(153) 

and 

1 ■ ■ 2 

9 ^ ^ = —A 9 ^ ( 29 ij 9 ]^i 9ii9jk 9ik9ji) ~ ~^-^Qki • 

(154) 

Hence 

C\kl = Pki - ‘^Qkl + 3 9 kl • 

(155) 

Consequently, the second invariant of Ting reads 



ft = c-= ft,P« - 4ft,Q« + 4Q„Q“ + 1 (S + 2 Af . 

(156) 

or 

B 2 = -Zl + -Z 2 + -Zs + Z 4 ^ — 4 Z 5 + 4:Zq . 

(157) 


A.3 The first Ahmad invariant 

For the first invariant of Ahmad, we need the trace 

C^kji = g^^Cikji = (^Sikji + • (158) 

Here, 

g^^Sikji = Ski = Pki + -^Sgki , (159) 

g ^ ^ ^Aikji = ~^Ag ^ {“^gikgji gugjk gijgki) — ~^Agki . (160) 

and 

g ^ ^ ^Aikji — g ^ {gijQki T gjkQu T guQjk T guiQij “^gjiQik ‘^gikQji) — Qki ; (161) 

Chji = Pki + Qki + -^gki{S - A). 


Consequently, 


( 162 ) 
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Hence using (155) and (162) we get 


Bs = {Pki + Qki + \9ki (5 - ^ (-5 + 2^) 


= PkiP^^ - PkiQ^‘ - 2QkiQ^^ + ^{S + 2A){S - , 


-ifcZ 


kl 


1 


or 


Bs — -^1 + -Z2 — -Z3 + Z4 — Z5 — 2 Zq . 


(163) 

(164) 


A.4 The second Ahmad invariant 


This invariant is obtained by the use of the formula (162), 


B 3 = {p^^ + + ^/'(5 - 41)^ Tpfcz + Qki + \gki{S - 4) 


1 


or 


= P«P“+ 2 P.< 3 »+Q.< 3 » + ^(S-A)^ 


-B 4 — -Zi — -Z 2 + 2^3 + -^4 + 2^5 + Zq 


(165) 

(166) 


A.5 The Norris invariant 

We put the invariant of Norris first in the form 

B^ = + AA') {S,kji + Akji) = 

Using ( 144|145|146), we have 


1 

cijkl c Qijkl Q ^ '7 \ ^ 'V \ '7 

* ^ikjl = * ^ijkl = + ^7 • 

We observe 

Thus we find 

^^Ajki = AA^{2gi,gki - gugjk - g^kg,l) {2g^^g^' - 

^ '^Aijki ^ ^4 ^ — (^gikQji T gjkQii T guQjk T gjiQik 2gkiQij 2gijQki) x 
[g"A^^ + g^^Q"^ + + g^A"^ - 2 gAQ^^ - 2 g^^QA) 

= -GQijQ^P 

B 5 = C^^^Aikji = ^Zi - ^Zs + ^Z4 - 6 Ze + ^7 • 


(167) 


(168) 


(169) 


(170) 


(171) 


Consequently, 


( 172 ) 
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